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Abstract 

Compactons are compactly supported solitary waves for nondissipative evolution 
equations with nonlinear dispersion. In applications, these model equations are ac- 
companied by dissipative terms which can be treated as small perturbations. We 
apply the method of adiabatic perturbations to compactons governed by the K (n, n) 
Rosenau-Hyman equation in the presence of dissipative terms preserving the "mass" 
of the compactons. The evolution equations for both the velocity and the amplitude 
of the compactons are determined for some linear and nonlinear dissipative terms: 
second-, fourth-, and sixth-order in the former case, and second- and fourth-order 
in the latter one. The numerical validation of the method is presented for a fourth- 
order, linear, dissipative perturbation which corresponds to a singular perturbation 
term. 

Key words: Soliton perturbation theory, Adiabatic perturbation method, 
Compactons, Solitons, K(n, n) Rosenau-Hyman equation 
AMS Codes: 35Q51; 35Q53; 35B20. 
PACS Codes: 02.30.Jr; 02.30.Mv. 



1 Introduction 



Generalized Korteweg-de Vries equations with nonlinear dispersion can propa- 
gate compactly supported solitary waves, referred to as compactons jl ,2,3 E1f5ll6lf 
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Numerical simulations show that an initial pulse wider than a compacton de- 
composes into a set of compactons with a small amount of radiation; more- 
over, compactons collide elastically suffering only a phase shift after the colli- 
sion and generating a small- amplitude, zero- mass, compact ripple [Tfgf5][TUfTTj . 
First discovered in the (focusing) K(n, n) Rosenau-Hyman equation for the 
modelling of pattern formation in liquid drops, compactons have several ap- 
plications in physics and science [Tj, such the pattern formation on liquid 
surfaces [12J, nonlinear excitations in Bose-Einstein condensates [13], the lu- 
brication approximation in thin films [14J, or even the pulse propagation in 
ventricle-aorta system [13]. The K(n,n) equation is also the continuous limit 
of the discrete equations of a nonlinear lattice [TfT6] and has been general- 
ized to higher dimensions [17] . Let us also note that the K{m,n) equation 
has also other kind of solutions, such as elliptic compactons [T|18fl9] . loop 
solutions [2D], dark solitons [21], peakons [22] and other solitary waves [2321] • 
Finally, let us remark that several generalizations of the K(m, n) equation 
have also been considered in the literature, for example, the inclusion of time- 
dependent damping and dispersion [23], or the addition of fifth-order disper- 




In applications, one- dimensional nonlinear evolution equations with solitary 
waves are usually obtained as the leading order term of a perturbative or 
asymptotic expansion of the solution of a more complicated mathematical 
model [271128] . The perturbation is adiabatic when the "mass" associated to 



adiabatic perturbation methods are referred to as soliton perturbation the- 
ory [33f 34|I35| . Solitary waves are usually robust to adiabatic perturbations, 
even if they are singular, when their effect is weak; in such a case, the pa- 
rameters of solitons and compactons, namely the velocity and the amplitude, 
slowly change in time. Let us note that the K(n, n) equation does not have a 
Lagrangian, hence the robustness of a compacton under perturbations cannot 
be studied by using the Lyapunov stability method developed by Dey and 
Khare [36] for the Cooper-Shepard-Sodano equation [2ll3] . 

Adiabatic perturbations methods have been applied to the compactons of the 
K(2, 2) equation with second- and fourth-order linear dissipative perturbations 
in Ref. [37J and numerical simulations show their good accuracy in determining 
the evolution of the velocity and amplitude of the perturbed compactons |38j . 
The damping of the amplitude of the compacton is accompanied by the gener- 
ation of trailing tails resulting from the conservation of its "mass" [39l40|41j . 
This paper extends the adiabatic perturbation theory to compactons of the 
K(n, n) equation under weak dissipative perturbations. Next section recalls 
the main properties of the perturbed K(n, n) equation. Section [3] applies the 
adiabatic perturbation method to five dissipative perturbation terms. The 
validity of the method for singular perturbations is illustrated by using nu- 
merical solutions in Section HI Finally, the last section is devoted to the main 




sion [26 . 



the solitary wave is conserved 




; in the context of solitons, such 
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conclusions. 



2 The perturbed K(n, n) equation 

Let us consider the K(m, n) equation given by 

« f + (A + M» = o, (i) 



where u(x,t) is the wave amplitude, x is the spatial coordinate, t is time, and 
the subindex indicate differentiation. This equation has at least two invariants 
given by 

M = h= f udx, and P = J 2 = / u n+1 dx, (2) 



referred to as "mass" and "momentum" in analogy with the K(m, 1) equation 
which has another invariant given by 

2 m+ 1 J 



referred to as "energy" [42J. Some authors refer to P as energy in the case of 
n = m, but in our opinion momentum is more appropriate. 

The perturbed K(n, n) equation is given by 

u t + {u n ) x + {u n ) xxx = eV{u), (3) 



where |e| 1 is a small parameter and V(u) is a function of u and its 
spatial and temporal derivatives. Under weak perturbations, the robustness 
of compactons results in the preservation of their shape only with a slow 
time evolution of their parameters; in such a case, by introducing a slow time 
r = et, the solution of Eq. ([3]) can be written as 

u(x, t, t) = u c (x, t, t) + Ut(x, t, r), (4) 



where u c corresponds to a compacton whose parameters evolve slowly in time 
and ut is the non-compacton part of the solution, usually a trailing tail. The 
evolution of the parameters of the compacton u c (x, t, r) of Eq. Q corresponds 
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to the ansatz given by 




l/(n-l) 



u c (x,t } r) = \ 2 7 ^ cos' , (r,) 



for |£| < nir/(n— 1), and u c (x,t,r) = otherwise, where £ = £(x,t,r) = 
x — c(r) t; hence the compacton's amplitude is {2nc(r)/(n + l)} 1 ^™ -1 ), and 
its velocity is c(r). Let us remark that the solution (JHJ) of Eq. (j3J) with £ = 
was first obtained for n = 2 and 3 in Ref. [I], and generalized in Ref. [3] 
for n G (1,3); the last condition is required for Eq. (J5]) to be a solution in 
the classical sense. Let us also note that the derivation of both invariants in 
Eq. d2J) is independent of the number of continuous derivatives at both edges 
of the compacton (jSJ) for n G (1, 3). 

For perturbations which do not preserve the mass of the compacton, the evo- 
lution in the slow time of the mass can be used to determine the evolution of 
c(t). The integration in space of Eq. fl3]) yields 



— J udx = e J V[u)dx. (6) 

-oo — oo 



dr 



Taking eV(u) = —e u, with |e | "C 1, and introducing u c given by Eq. (JSJ) 
into Eq. yields 

_ = -(n-l)c(r), 



whose solution is 

c(r) = c (0)e- £o(n - 1)r , (7) 



indicating that both the velocity and the amplitude of the compacton decrease 
slowly as it evolves in time for n > 1. 



3 Adiabatic perturbation theory for compactons 



The adiabatic perturbation theory is a technique for the analysis of solitary 
wave solutions of nonlinear evolution equations under dissipative perturba- 
tions preserving the mass of the compacton (let's notice that in such a case 
the right-hand side of Eq. ([6]) is null). This technique determines the slow 
time evolution of the parameters of the compacton by means of using the slow 
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temporal derivative of the invariant of Eq. fl3]) with e = 0. Multiplying Eq. 
by u n and integrating in space results in 



oc 



d r u n+1 
dr J n + 1 



oc 

dx = J u n V{u) dx, 



where the substitution of u by u c in Eq. (jSJ) gives an ordinary differential 
equation for the velocity c(r) as function of the slow time r. In such a case, 
the integral of the left-hand side of Eq. (JSJ yields 



/ — n 2 " P I 3n+l | 

(n + i) U + i/ r(a±i) ' 

and its right-hand side results in a similar expression depending of the per- 
turbation V(u) considered. 

Next subsections illustrate this method by using several dissipative perturba- 
tions, both linear and nonlinear. Other perturbations not considered in this 
paper may be dealt with similarly. 

3.1 Perturbation with second-order derivatives 

Let us study the effect of a linear perturbation given by 

eV(u) = a u xx - ai u xt + a 2 u t t, (9) 



u, 



\n+l 



n + 1 



dx 



where < a, <C 1. By means of straightforward integration using standard 
mathematical software, Eq. (jSJ gives 

c'(r) = - 2 ^~^ffi^ (a c(r) + a, c(r) 2 + a 2 c(r) 3 ) . (10) 
n (n + 3) r + 2) v ' 

By using the properties of the Gamma function, for 1 < n < 3, this equation 
can be simplified to yield 



/( s (n-1) 



2 



n (n + 3) 



(a c(r) + a 1 c(r) 2 + a 2 c(r) 3 ). (11) 



The analytical solution of this equation can be obtained in implicit form, for 
a 2 = also in explicit form, showing that the perturbation <Q is dissipative 
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for aj > 0. 



3.2 Perturbation with fourth- order derivatives 



Let us take a fourth-order linear perturbation given by 



eV{u) = -P u 



xxxx 



l^xxxt - 02 Uxxtt + @3 Uxttt ~ Pi Utttt 



(12) 



where < -C 1. After a long integration process of Eq. (jHJ) gives 



(n - l) 3 {{n - 3) n - 1) 
{n — 5) n 3 (n + 3) 



(/3 c(r)+A c(r) 2 + (3 2 c(r) 3 
+ /3 3 c(r) 4 + /3 4 c(r) 5 ), 



for 1 < n < 3. The analytical solution of this equation can be obtained in 
implicit form and their numerical evaluation shows that the perturbation ffl2l) 
is dissipative for > 0. 



3.3 Perturbation with sixth-order derivatives 

The last linear perturbation to be considered in this paper is 

E 7^[U) 70 Hxxxxxx 7l ^xxxxxt 72 ^xxxxtt 73 ^xxxttt 
+ 74 Mraitft - 75 Mxtttit + 76 Utttttt, 

where < ji -C 1. For this perturbation, Eq. (jHJ) yields 



for 1 < n < 7/3; note that the integration of Eq. (jHJ) for n > 7/3 does 
not converge to a finite value. The analytical solution of this equation can 
be obtained in implicit form and their numerical evaluation shows that the 
perturbation ( fT2l) is dissipative for 7, > 0. 




+ 72 c(r) 3 + 73 c(r) 4 + 74 c(r) 5 + 7s c(r) 6 + 7e c(r) 7 ) , 
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3.4 Nonlinear perturbation with second- order derivatives 



Dissipative perturbations can also be nonlinear, for example, the second-order 
nonlinear perturbation given by 

£ V(u) = 61 (n - 1) n {u x f u n ~ 2 + 5 2 n u xx u n ~\ (13) 



where < 5i <C 1. Note that eV{u) = 5 (u n ) xx for 8\ = 5 2 = 6. By means of 
straightforward integration, in this case, Eq. OH]) gives 

c'(r) = (5 1 (n - 1) + 8 2 (1 - 2 n)) ^ J c(r) 2 . (14) 



This equation reduces to 



for 5i = 5 2 = S. 

The analytical solution of Eq. ( |T5|) shows that the perturbation ( |T3|) is dissi- 
pative for 5 > 0, but that of Eq. ( 1T41) shows that a dissipative perturbation 
requires, for 52 = 0, that 5\ < for 1 < n < 2, or 5\ > for 2 < n < 3, 
and, for <5i = 0, that 5 2 > for 1 < n < 3/2 or 2 < n < 3, or <5 2 < for 
3/2 < n < 2. 

5.5 Nonlinear perturbation with fourth- order derivatives 
Our final example of a nonlinear dissipative perturbation is 

eV{u) = m {u x f u n - A + V2 {u x f u xx u n ~ 3 + m (u xx ) 2 u n ~ 2 

-\- 774 U x U xxx U -\- T]^ U xxxx U , 

where < rji <C 1. After a long integration, Eq. (jSJ) yields 



c'(r) 



(n - l) 3 c(r) 2 
2n 4 (n + T) (n + 3) 



3 r/i - 772 (2 n - 3) + r? 3 ( 2 n 2 - 2 n + 3 



+774 (2 n 2 - 8 n + 3) - r) 5 (2 n - 1) (2 n 2 - 8 n + 3)) . (16) 
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The conditions on rji such that the solution of this equation is dissipative must 
be discussed in a term- by-term basis with a r/j 7^ and the other ones rjj = 0, 
for j 7^ i. The first term is dissipative for 7]i < 0; the second one, for 7/ 2 < 
with 1 < n < 3/2, and for 772 > with 3/2 < n < 3; the third one, for 773 < 0; 
the fourth one, for 774 > 0; and the fifth one, for 775 < 0. Further possibilities 
require a straightforward but cumbersome combination of these facts. 

Note that for eV(u) = rj (u n ) xxxx , i.e. for 7]i = 77 (n — 3) (n — 2) (n — 1) n, 
V2 = 6 77 {n — 2) {n — 1) n, r] 3 = 3 r\ {n — 1) n, 774 = 4 7/ — 1) n, and r/ 5 = r/n, 
Eq. fTl6|) reduces to 



whose analytical solution it dissipative for r] < 0. 



4 Validation for singular perturbations 

The main criticism on the validity of the adiabatic perturbation method pre- 
sented in the last section is the question of its applicability to singular pertur- 
bations, i.e., when eV(u) has more than three spatial derivatives. Such terms 
may cause a fundamental change at the edge of the compacton and may in- 
validate the whole structure. Let us check by numerical simulation that, for 
small enough e, the method works properly. 

The numerical solution of the perturbed K(n, n) equation with a fourth-order 
linear, dissipative term given by 

u t -c u x + {u n ) x + (u n ) xxx + A) u xxxx = 0, (18) 



can be obtained by the Petrov-Galerkin method presented in Refs. [8|38|10j . 

using the implicit midpoint rule for the integration in time [H]; a detailed 
presentation is omitted here for the sake of brevity. 

The perturbed compacton @ develops a trailing tail under the dissipative 
term, but evolves preserving the mass of the solution, cf. Eq. (T5]), hence 

M = J u(x, t, r)d X = J u c ( X , t, r) dx + J M*, «, r) dx. (19) 
From this expression the area under the tail can be easily calculated yielding 



c'(r) = 7] 



(n- l) 3 (2 + n) c(r) 2 
2n(n + l) (n + 3) 



(17) 




8 




0.4 



K(2,2) 



3 







500 1000 1500 



2000 



Fig. 1. Evolution in time of the area under the tail of a compacton calculated by a 
quadrature rule applied to the numerical solution (solid line) and with an adiabatic 
approximation (dashed one) for the K(n,n) equation with n = 2, 5/3, 3/2, 7/5, 
4/3, 9/7, and 5/4. The compacton has a velocity c = 1, and the numerical method 
uses c = !,#,= 0.001, Ax = 0.2, and At = 0.1. 



Figure [T] shows the evolution in time of the area under the tail of a compacton 
calculated by a quadrature rule applied to the numerical solution of Eq. (JIB]) 
(solid line) and by using Eq. ( 120]) (dashed one) for the K(n, n) equation with 
n = 2, 5/3, 3/2, 7/5, 4/3, 9/7, and 5/4. The plot shows the good accuracy of 
the adiabatic perturbation method even for a singular perturbation. 

Equation f[2"0"]) can be used to estimate the shape of the trailing tail, except 
at its front. The left edge of the perturbed compacton is located at position 
X(r) given by 



where X(0) is the location of the left edge of the compacton at t — 0. Hence, 
the tail area (1201 can be calculated as 

X(t) 

At(t) = J UT(x,t,r) dx. 

X(0) 

The shape of the tail can be obtained by using Leibniz's rule for differentiation 




ut(x, t, r) dx = M 



(20) 
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Fig. 2. The shape of the tails of a compacton obtained by means of the adiabatic per- 
turbation method (solid line) and calculated by using a numerical method (dashed 
one) for the K(n,n) equation with fourth-order linear dissipations /?o = 0.01, 0.005, 
and 0.001, for n = 3 (left top plot), 2 (right top one), 7/5 (left bottom one), and 5/4 
(right bottom one). In all the plots the compacton has a velocity c = 1, t = 1900, 
and the numerical method uses cq = 1, Ax = 0.2, and At = 0.1. 

under the integral sign, yielding 

M X(r), f ,r) = ^^. (21) 



Figure [2] shows the shape of the tails of a perturbed compacton obtained by 
means of Eq. ff2T]) (solid line) and numerically calculated by using Matlabs's 
ODE suite (dashed one) for Eq. ffLSl) with n = 3 (left top plot), 2 (right top 
one), 7/5 (left bottom one), and 5/4 (right bottom one); each plot presents 
three curves for (3 = 0.01, 0.005, and 0.001. The plots zoom into the tail 
because of its small amplitude compared to that of the compacton. For the 
K(2, 2) (right top plot) the agreement between the adiabatic perturbation 
method and the numerical result is extremely good due to the small width of 
the front of the trailing tail. As the value of n < 2 approaches unity from above, 
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the width of the front of the tail increases, resulting in a loss of the accuracy 
of Eq. (l2Tj) as shown in the bottom plots in Fig. [2] even in such for 
small enough values of the parameter /3 , the agreement between perturbation 
and numerical results is good. For the limiting case n = 3 (left top plot in 
Fig. [2]), the accuracy of the adiabatic perturbation method developed in this 
paper is worst except for small values of /3q. 



5 Conclusions 

The adiabatic perturbation method has been applied to the K (n, n) Rosenau- 
Hyman equation with both linear and nonlinear dissipation terms. The slow 
time dynamics of the compacton velocity (which uniquely determines that 
of its amplitude) has been determined for five different perturbations. The 
analytical results have been validated by means of numerical methods for a 
singular perturbation, showing the good accuracy of the adiabatic perturba- 
tion method, even for the estimation of the shape of the trailing tails of the 
perturbed compactons. 

The adiabatic perturbation method can be applied to other nonlinear evolution 
equations under dissipative perturbations having compacton solutions, such 
the K*(n, n) Cooper-Shepard-Sodano equation, and generalized versions of the 
Boussinesq, regularized long- wave, Benjamin-Bona-Mahony, and Camassa- 
Holm equations, to mention only a few [3H] . 
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